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Abstract 

£SJ ' Narrow-line Seyfert 1 galaxies (NLSls) exhibit extreme soft X-ray excess and large variability. We 

argue that both features can be basically accounted for by the slim disk model. We assume that a central 
black-hole mass in NLS1 is relatively small, M ~ 10 5_7 Mq, and that a disk shines nearly at the Eddington 
luminosity, Lg- Then, the disk becomes a slim disk and exhibits the following distinctive signatures: (1) 
The disk luminosity (particularly of X-rays) is insensitive to mass-flow rates, M, since the generated energy 
is partly carried away to the black hole by trapped photons in accretion flow. (2) The spectra are multi-color 
blackbody. The maximum blackbody temperature is Xbb — O.2(M/1O 5 M0) -1 / 4 keV, and the size of the 
blackbody emitting region is small, rbb % 3rg (with r$ being Schwarzschild radius) even for a Schwarzschild 
black hole. (3) All the ASCA observation data of NLSls fall onto the region of Af/(L E /c 2 ) > 10 (with L E 
being the Eddington luminosity) on the (rt,b>2bb) plane, supporting our view that a slim disk emits soft 
X-rays at ~ Le in NLSls. (4) Magnetic energy can be amplified, at most, up to the equipartition value 
with the trapped radiation energy which greatly exceeds radiation energy emitted from the disk. Hence, 
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energy release by consecutive magnetic reconnection will give rise to substantial variability in soft X-ray 



emission. 



Introduction seems that a group of NLSls is smoothly connected to 
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H the class of BLSls (Brandt, Boiler 1998). 

The terminology of narrow-line Seyfert 1 galaxies It is often suggested that the strength of the soft X-ray 

(NLSls) originates from their relatively narrow optical excess relative to the hard X-ray power law appears to be 

Balmer-line emission with v B lr ^ 2000km s" 1 compared directly related to the 'primary eigenvector' of Boroson 

with those of usual Seyfert 1 galaxies with broad-line and Green (1992), which represents the strongest set of 

emission, BLSls (see a concise review by Brandt 1999 optical emission line correlations (e.g., Brandt 1999). In 

and references therein). Their optical line ratios and re- this sense ; NLSls are located at the extreme end of the 

markable Fen lines (Osterbrock, Pogge 1985; Halpern, primary eigenvector. The UV line properties of NLSls 

Oke 1987; Grupe et al. 1999) are also distinct from those also fits this scheme (Wills et al. 1999). Then, what is the 

of BLSls. NLSls exhibit unique X-ray properties {e.g., control parameter which drives the primary eigenvector? 

Grupe et al. 1998); they are characterized by large soft What physical factor distinguishes NLSls from BLSls? 

X-ray excess (Pounds et al. 1996; Otani et al. 1996; There are several models proposed to account for the 

Leighly 1999b) and a good correlation is known to exist narrowness of otherwise broad-line emission (see Boiler 

between the strength of soft excess and FWHM of optical et al. 1996 for extensive discussion). Successful models 

Balmer lines (Boiler et al. 1996; Laor et al. 1997). Rapid should account for a smooth continuation between NLSls 

soft/hard X-ray variability is another signature charac- and BLSls. Among them, one fascinating and probably 

terizing NLSls (Otani et al. 1996; Boiler et al. 1997; most promising explanation is that NLSls contain rel- 

Leighly 1999a). NLSls are not rare but comprise a sig- atively less massive black holes (with M ~ 1O 5-7 M0). 

nificant part (say, 20%) of the Seyfert 1 galaxies and it Moderate brightness of NLSls, L ~ 10 43_45 erg s _1 , im- 
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plies similar mass-flow rates, M, to those of BLSfs. It 
then follows that the ratio of L/L-e (or M/M) is rel- 
atively large for NLSls. The control parameter which 
drives the eigen vector could be the fraction of the Ed- 
dington rate at which the supermassive black hole is ac- 
creting, M/(L E /c 2 ). 

For small black-hole masses, narrow Balmer line emis- 
sion of NLSls can be understood, provided that the 
broad-line clouds (which emit broad-line emission in 
BLSls) are bound in the potential by a central black 
hole (Laor et al. 1997). Since the radius of the broad- 
line clouds are roughly scaled as tblr oc L - 5 from the 
reverberation mapping (Kaspi et al. 1996), the circular 
velocity of broad- line clouds, ~ GMj Vblr, is system- 
atically smaller for smaller M at a constant L (and thus 
constant tblr), yielding narrow Balmer emission. 

For a luminosity close to the Eddington, 

L E ^ 1.2 x 10 43 M 5 erg s" 1 , (1) 

where M 5 = M/(1O 5 M ), what do we then expect the- 
oretically for the disk structure and its spectra? It is 
known that for such a high luminosity advective energy 
transport dominates over radiative cooling (Abramowicz 
et al. 1988; section 10.3 of Kato et al. 1998). Such 
a disk is named as the slim disk, since it is moderately 
geometrically thick. 

The slim disk is an optically thick version of ADAF 
(advection-dominated accretion flow) and should not be 
confused with optically thin ADAF proposed for low- 
luminosity AGNs (see section 10.2 of Kato et al. 1998 
and references therein). Unlike the optically thin ADAF 
the observable features of the slim disk have been poorly 
investigated. Szuszkiewicz et al. (1996) was the first to 
start investigation along this line (see also Wang et al. 
1999), and we give more detailed discussion on the ob- 
servational consequences in the present study. Watarai 
et al. (2000a) studied the super-critical accretion flow in 
the context of Galactic black-hole candidates (GBHC), 
showing some unique observational features of the slim 
disk. 

We focus on the AGN case in the present study. Meth- 
ods of calculations are given in section 2. We then discuss 
the non-standard emission properties of the slim disk in 
the context of NLSls in section 3. We next turn to the 
subject of large variability and discuss its origin in rela- 
tion to magnetic-field activity in section 4. Outstanding 
issues will be discussed in section 5. The final section is 
devoted to conclusions. Throughout the present study 
we use the normalization of M with M cr n = Le/c 2 ; i.e., 

m = M/M crit ~ M/(1.3 x 10 22 M 5 g s" 1 ) 

- M/(2.5 x 1Q- 4 M 5 M G yv- 1 ). (2) 



2. Basic equations and methods of calculations 

We first solve the steady-state, transonic disk struc- 
ture. The methods are the same as those adopted in 
Matsumoto et al. (1984) and Watarai et al. (2000a), but 
we briefly repeat the description below. 

For the purpose of the present study, it is essential 
to solve the full-set equations, since the emission from 
the innermost region is in question. (Approximate so- 
lutions, such as the self-similar solutions, do not prop- 
erly treat the inner boundary conditions nor regular- 
ity conditions.) The transonic nature of the accretion 
flow should also be carefully treated. We do not include 
any general relativistic (GR) effects except for adopting 
the pseudo-Newtonian potential by Paczynski and Witta 
(1980), V = ~GM/{r-r s ). For discussion of GR effects, 
see section 3.5. We use cylindrical coordinates, (r, ip, z). 

Throughout this paper, we used vertically integrated 
(or height-averaged) equations. In fact, we numerically 
confirmed that H/r < 1 for rh< 100, where H is the 
scale-height of the disk. Certainly, H/r > 1 for even 
higher mass accretion rate, but for such high rh, pho- 
ton diffusion time over H becomes longer than accretion 
timescale so that the disk height should be significantly 
reduced (see section 5.1). The height-average approach 
may not be very precise and we eventually need fully 
2-D models. The same problem arises in optically thin 
ADAF, for which Narayan & Yi (1995) found that the 
2-D solutions of the exact non-height-integrated equa- 
tions agree quite well with those of the simplified height- 
integrated equations, provided that 'height-integration' 
is done along 6 at constant spherical radius, rather than 
along z at constant cylindrical radius. The height- 
integrated equations therefore are a fairly accurate rep- 
resentation of quasi-spherical ADAFs. 

We construct the vertically integrated equations, using 
integrated variables, such as surface density, S = J pdz, 
and integrated pressure, II = J pdz. Vertically integrated 
mass, momentum, and angular-momentum conservations 
lead 

- 27rr£u r = M = const., (3) 
dv r 1 dU _ £ 2 - i\ II dlnft K 

and 

M{£ - £ in ) = -2irr 2 T rv , (5) 
respectively, 

where fi(= v v /r) and fi K [= (GM/r) 1 / 2 /^ - r s )\ are 
the angular frequency of the gas flow and the Keplerian 
angular frequency in the pseudo-Newtonian potential, £ 
(= rv v ) is specific angular momentum, (= r 2 rix) is 
the Keplerian angular momentum, £- ln is the specific an- 
gular momentum finally swallowed by the black hole, and 
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r s = 2GM/c 2 ~ 3 x 10 10 M 5 cm is the Schwarzschild 
radius. As to the viscous stress tensor we adopt the 
prescription T rtp = —all with the viscosity parameter, 
a = 0.1. Note, however, that the resultant disk structure 
does not sensitively depend on a (Watarai et al. 2000b). 

The vertically integrated energy equation is symboli- 
cally written as 



Qadv ~ Qvis Qrad' 



(6) 



where Q~ dv [oc Hv r T(ds/dr)] is the advective cooling with 
s being specific entropy, and other two terms on the right- 
hand side represent viscous heating and radiative cooling, 
respectively. 

The outer boundary conditions are imposed at r = 
2.0 x 10 3 rs, where each physical quantity is taken as that 
of the standard disk. We confirmed that the disk struc- 
ture near the outer boundary can be very accurately de- 
scribed by the standard-disk relation. When calculating 
spectra (see below) , we extrapolate the temperature dis- 
tribution up to r — 2.0 x 10 4 rs using the standard-disk 
relation, T e gf cx r~ 3 / 4 . The basic equations are integrated 
by the semi-implicit method from the outer boundary to 
the inner one taken at r = l.Olrs. The flow is subsonic 
far outside (r ^> rs) but is supersonic just near the black 
hole. The solution should satisfy the regularity condition 
at the transonic radius, at about 2.7rg, and we adopt the 
free boundary conditions at the inner edge. 

3. Observational features of a slim disk 

3.1. Unique temperature profile 

How do structural changes occur from the standard 
to the slim disk when rh increases? The solid lines in 
the upper panels of figure 1 represent the temperature 
profiles of the disks for various values of rh for the black- 
hole mass of M — 10 5 M Q . For comparison, we also plot 
with the dashed lines the temperatures expected from 
the standard-disk relation (Shakura, Sunyacv 1973), 



T eff ~ 10 6 5 M, 



rs 



1 




(7) 



where M 5 = M/10 5 M Q . 

We note three characteristics: 

1. In contrast with the standard disk having an effec- 
tive temperature profile of T c s oc r -3 / 4 , the slim 
disk exhibits flatter slopes, T e g oc r -1 / 2 (Watarai, 
Fukue 1999). As rh increases, the part with a flat- 
ter temperature profile expands outward; i.e., up 
to r ~ 100rs for rh ~ 100, and r ~ 1000rs for 
m - 1000. 

2. When the innermost parts become advection- 
dominated at rh > 50, substantial amount of radi- 
ation is expected from inside the marginally stable 
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Fig. 1.. The temperature profiles (upper) and the emergent 
spectra (lower) of the calculated disks with different mass-flow 
rate; m = 1.0, 10.0, 100.0, and 1000.0 from the below. As 
rh increases, regions with a fiat temperature profile (extending 
down to logr ~ 10.5 or r ~ rg) expands from the innermost 
region and the disk spectra become flatter at the peak. For 
comparison, we plot the temperature profiles expected by the 
standard-disk relation [equation (jij)] and their emergent spectra 
by the dashed lines. 



last stable orbit, r ms = 3rg ~ 10 11 cm. This can be 
understood as follows: For rh -C 1, the flow is very 
subsonic far outside, while it should become super- 
sonic near the black hole. A sudden jump in the 
radial accretion velocity across the transonic radius 
(at ~ 2.7rs) results in an abrupt density drop inside 
the transonic radius (note M oc u r E is kept con- 
stant) . The innermost region become of low density, 
thus producing little emission (emissivity goes like 
p 2 ). For m > 50, conversely, the flow velocity out- 
side the transonic radius is already comparable to 
the sound speed, thus the velocity and density jump 
across the transonic radius being small. Therefore, 
huge, continuous mass supply makes it possible to 
keep the innermost region (at r < 3rs) optically 
thick and substantial amount of radiation can be 
produced there as a result. 

3. Consequently, the inner disk temperature (which 
makes a hotter peak in the spectrum) gets system- 
atically higher than the values expected by the stan- 
dard relation. The size of the region emitting with 
the maximum energy becomes small, however. 
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3. 2. Spectral features 

Once temperature distributions are obtained, we can 
calculate emergent spectra by simply summing up a con- 
tribution from each tiny portion of the disk. Such scatter- 
ing dominated, hot disks as those considered here do not 
always emit blackbody radiation (Czerny, Elvis 1987; Be- 
loborodov 1998). However, we simply assume blackbody 
radiation to see the effects of transonic flow (discussed 
later). 

The results are plotted in the lower panels of figure 1. 
As to increases, the spectra, vS v , get flatter and flatter 
around the peak, vS v ~ v°. (Still, vS v ~ J/ 4 / 3 in the 
optical bands.) According to Grupe et al. (1998) soft X- 
ray selected AGNs (including a number of NLSls) tend to 
exhibit bluer optical/UV (around 1400 A) spectra when 
compared with those of hard X-ray selected ones and the 
maximum slope is S v oc v 3 . This is consistent with our 
calculations. 

The most significant problem lies in a os , or the opti- 
cal to soft X-ray flux ratio (say 3000A to 0.2 keV). High 
L/L-e models predict a very flat a os , -1 or flatter, as 
shown in figure 1 , whereas the observed slope is more like 
-1.5 (see Table 1 of Brandt et al. 2000, which includes 
many NLSls). Thus, the observed soft X-ray emission is 
a factor of 10, or more, weaker than high L/L-e disk mod- 
els predict. In fact, this problem already exists for normal 
AGNs (Laor et al. 1997) and is actually more severe for 
NLSls because of the flatter predicted a os . This does not 
necessarily rule out the disk model, but it does suggest a 
significantly reduced efficiency of the emission from the 
inner parts of the accretion disk (e.g. Laor 1998). 

NLSls have a steeper soft (0.1-2 keV) X-ray slope, but 
since they have a pretty normal a os , this slope is actually 
mostly because the hard X-ray component is weak and 
not because the soft component is strong (Laor et al. 
1997). The slim disk model shows that exponential roll- 
over extends over IkeV, thus producing large power-law 
indices in the ROSAT band. The inclusion of Compton 
scattering within the disk will increase the photon indices 
further (Ross et al. 1992) which will be discussed in 
section 3.5. Fairly strong correlation between the optical 
and soft X-ray emission is naturally explained since both 
are from the same optically thick disk. 

If we include a hard power-law component in the cal- 
culation, we expect a larger hard power law index in the 
ASCA bands for a larger soft excess as a result of more 
efficient cooling of electrons via inverse Compton scat- 
tering. Since the hard power-law slope is determined by 
the Compton y-parameter which is related to the energy 
amplifications of soft photons (e.g. Rybicki, Lightman 
1979), less (or more) energy input to high-energy elec- 
trons compared with that to soft photons leads to a larger 
(smaller) y and thus a larger (smaller) photon index in 
the hard X-ray bands; a x ft 2 (a x < 2) as is observed in 
the soft (hard) state of GBHC (Mineshige et al. 1995). 
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Fig. 2.. Disk luminosity as a function of rh. Asterisks denote 
the calculated luminosities, whereas the solid line shows the 
fitting formula [equation (H)]. It is clear that increase in L is 
suppressed at L > 2Lg. 



Systematically large a x 's in NLSls in the ASCA bands 
(Pounds et al. 1996; Brandt et al. 1997) can be under- 
stood. 

Further, Tanaka (1989) found highly variable hard 
emission relative to the soft one in GBHC. If the analogy 
holds, the relative strength of the hard component could 
take any values. This may explain a diversity in a x 's of 
NLSls in the ASCA bands. 

3.3. Evolution of the disk luminosity 

It should be repeated that the disk luminosity, in- 
creases only moderately above Le for steadily increas- 
ing M. We calculate L by summing up the contribution 
from each concentric annulus and plot in figure 2 as a 
function of rh. Asterisks are calculated values, while the 
solid curve represents the fitting formula (see Watarai et 
al. 2000a for a derivation), 



L(rh)/Li 



2[l + ln(m/50)] for rh > 50 
to/25 for rh < 50 



(8) 



The functional dependence of L /Le on rh is basically the 
same as that for GBHC with smaller black-hole masses, 
thus being insensitive to M. 

Care must be taken when interpreting L, since ra- 
diation field is likely to be anisotropic; i.e., radiation 
mostly goes out in the vertical direction to the disk plane, 
whereas matter approaches to the black hole along the 
disk plane. This is a way in which super-critical accretion 
is possible, but at the same time, we expect much diluted 
radiation going out in the direction of the disk plane. 
Blocking of the radiation by the outer disk rim provides 
an additional diluting effect (Fukue 2000). Thus, we 
will detect weaker flux (< Le/^oI 2 with d being dis- 
tance), even when the total L exceeds Le- We also ex- 
pect strong radiation-accelerated wind in the central re- 
gion (e.g., Watarai, Fukue 1999), which will give rise to 
blue-shifted UV lines (Leighly 1999b). 
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3.4- Fitting with the multi-color blackbody 

Now we fit the calculated disk spectra by the extended 
disk blackbody model (Mineshige et al. 1994a), assuming 
the following T c g profile, 

T cff = T in (r/r in yP, (9) 

where r; n , T; n , and p are fitting parameters. 

Note that the fitting results may slightly depend on the 
X-ray energy ranges for which fitting is made, but in the 
present study we fix the energy bands to be either 0.1 - 
0.5 keV (for m < 10), 0.1 - 1.0 keV (for 10 < m < 100), 
or 0.1 - 3.0 keV (for rh > 100). 

To have a good fitting result, especially in p, wider 
spectral ranges, especially a good coverage of soft X-ray 
bands, are certainly necessary. The small black hole mass 
is preferable for this sort of fitting, since then the disk 
temperature is roughly scaled as <~ M -1 / 4 , expected from 
the standard-disk relation (which is also roughly valid 
in the slim-disk regime). In the fitting, T; n and n n are 
basically determined by the energy and the intensity of 
the exponential roll-over, while p depends on the slope in 
lower energy bands. EUV to soft X-ray observations are 
crucial for discriminating the slim disk from the standard 
disk. 

We examine how the derived fitting parameters of r- m , 
T; n , and p change with an increase of a mass-flow rate 
(and L). The results of the fittings are summarized in 
table 1 and figure 3. 

As m increases, so does T; n more or less gradually, 
but at the same time n n shifts from <~ 10rs to ~ 0.5rs 
abruptly at L/Te ~ 2 (m ~ 50). We should recall that 
r; n stands for the size of the area emitting with B v {T- lri ) 
and not the radius of the physical inner boundary (which 
should be always at rs). In fact, the upper part of figure 
1 has already shown that the peak effective temperature 
shifts toward smaller radii as rh increases. 

The effective temperature of the slim disk has a weak 
luminosity dependence, T; n oc M -1 / 4 cx L -1 / 4 , however, 
Tj n also has some dependence on M/M at rh < 50. We, 
hence, expect scatters in Tj n even for the same M and L. 



Table 1. Results of fitting. 
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Fig. 3.. The results of the fitting, r- ln (by the solid lines) and Ti n 
(by the dotted lines) as functions of rh (upper) and of L/L^, 
respectively. 

3.5. Comparison with the observations 

Hayashida (2000) analyzed the ASCA spectral data of 
9 NLSls and fitted them with a two spectral-component 
model: a hard power-law component and a soft black- 
body one. We assume that the soft component repre- 
sents the emission from the disk we are discussing. The 
fitting parameters of the soft component are blackbody 
temperature, Tbb, and the size of the blackbody emitting 
region, rbb. The results of the fitting are summarized in 
figure 4 and in table 2. 

Note that the redshift are small, z < 0.1, for all the 
sources. Among the nine sources in Hayashida (2000), 
we excluded Mrk766 because of its complex spectral fea- 
tures and their variability (Lcighly et al. 1996). Even 
for the other eight sources, it must be kept in mind that 
there may be some model dependence in the data points 
of Tbb and in rbb- In particular, modeling of peculiar 
spectral features around 1.1 keV found for some sources 
[IRAS13224, H0707, and PG1404; see Hayashida (2000) 
and reference therein] affect the parameter Tbb signifi- 
cantly, e.g., T bb for H0707 might be about 30% higher 
than the value listed in table 2, if we model the feature 
as absorption edges. On the other hand, sources of which 
soft component arc not so prominent, e.g., IZwl, the re- 
sults will be affected, if we allow some deviation from the 
perfect power law model for the hard component. Cal- 
ibration problem of ASCA SIS may add extra error on 
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the results, too. These points would be addressed in a 
separate paper. 

The theoretical expectations based on the slim disk 
model are also plotted by the solid and dashed lines in 
figure 4. Since inclination angles are not known, we as- 
sume the inclination angle of i = 60° in the theoretical 
models; i.e., we set r bb = r in Vcosi. Likewise, we set 
Tbb = 7i n , although the former may be systematically 
higher than the latter, T; n ~ (1.1 — 1.3)Tbb- Accordingly, 
the value of rbb should be smaller, for L oc r^ h T^ h is 
kept roughly constant. It might also be noted that the 
combination of rbb^bb' which is proportional to L/Le, is 
mass-independent, since rbb (oc r; n ) can be scaled with 
rs oc M and Tbb ~ T n oc M -1 / 4 . Therefore, constant m 
lines give a relation, Tbb oc r^^ 4 . 

On the theoretical grounds, various kinds of GR effects 
and Compton scattering should affect the results. When 
scattering opacity dominates over absorption, the emer- 
gent spectra deviate from blackbody one (see Rybicki, 
Lightman 1979). In a hot disk with Compton y > 1, 
moreover, inverse Compton scattering should modify the 
disk spectra in the high-energy tail, thereby yielding, at 
most, a factor of 2 increase in Tbb, while the total flux 
(oc TbbTbb) is kept constant; i.e., a factor of 4 decrease 
in rbb (Ross et al. 1992; Shimura, Takahara 1995). 

GR effects (such as gravitational redshift, Doppler 
boosting, gravitational focusing due to ray bending) pro- 
duce complex effects. We calculate a slim disk model 
with face-on geometry, for which only gravitational red- 
shift works, finding systematically larger rbb (~ 3r§) by 
a factor of ~ 6 and thus lower Tbb (~ 0.12 keV) by a fac- 
tor of ~ 2.5 even for m > 100. In other words, most of ra- 
diation from inside 3r$ vanishes for this case. [Although 
photon energy gets lost partly, the total flux does not 
drop much because of flatter T e g profile, T e ff(r) oc r -1 / 2 .] 
Since Compton and GR effects work independently in op- 
posite ways for a face-on disk, these effects cause decrease 
in Tbb (by a factor of 2.5/2) and increase in rbb- Correc- 
tions to remove Compton and GR effects from the data 
points are, hence, A log Tbb < +0.1 and A log rbb < —0.2, 
which is indicated by the lower-right- ward arrow in figure 
4. (Since a model predicts lower Tbb by these effects, we 
need higher m to reproduce the observed Tbb-) 

For disks with non-zero inclination angles around a 
Schwarzschild black hole, Doppler boosting enhances ra- 
diation, yielding, at most, a factor of ~ 1.4 increase in 
Tbb compared with face-on disks, while the total flux kept 
roughly constant (Sun, Malkan 1989). We thus expect 
slight increase in Tbb and decrease in rbb for a fixed 
r bb-^bb compared with face-on disks. Required correc- 
tions are A log Tbb ^ — 0.2 and A log rbb ~ +0.4, which 
is indicated by the upper-left-ward arrow in figure 4. 

It is evident in figure 4 that all NLSls have relatively 
large m. In other words, NLSls shine at luminosities 
close to the Eddington (Hayashida 2000). Especially, 
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Fig. 4.. The r bb -T bb diagram of the observational data of NLSls. 
The loci of fa = 10 and 100 (by the solid lines) and those of 
m(= M/Mq) = 10 5 and 1O 6 M (by the dotted lines) obtained 
by the slim disk arc also shown (i = 60° is assumed). The arrow 
headed for upper-left (or lower-right) from the filled circle in 
the lower-left corner indicates the direction of the correction 
for that point to remove Compton and GR effects in the case 
of a face-on (nearly edge-on) disk. 

PG1404, H0707, and IRAS13224, all having large m in 
figure 4 and enhanced soft components, are known to 
exhibit giant amplitude variability. In fact, Hayashida 
(2000) found that the mass estimated by the variabil- 
ity (Hayashida et al. 1998) significantly deviate from 
that estimated by equating rbb = 0.5rs in those sources, 
suggesting anomalously large variability amplitudes com- 
pared with other NLSls with a similar black-hole mass. 
This tendency indicates a correlation between L/Le and 
the variability amplitudes (and soft X-ray excess over 
hard X-rays, see below). 

Despite ambiguities in data points in figure 4, we may 
conclude that all (or, at least, most of) the NLSls fall on 
the region with m > 10, supporting our view that slim 
disks emit substantial fraction of soft X-rays at ~ Le 
in NLSls. The observed steep slopes, a os (see section 
3.2), rather enhances our conclusion of high m, since it 
means more optical radiation and thus higher m than 
that expected solely from X-rays. 



Table 2. Results of ASCA Observations. 



source/m 


T bb (keV) 


r bb (cm) 


log(T b 4 b r bb ) 


IZwl 


0.100 


9.75 x 10 1() 


6.99 


REJ1034 


0.107 


1.04 x 10 11 


7.14 


PG1440 


0.097 


2.08 x 10 11 


7.28 


IRAS13224 


0.098 


2.55 x 10 11 


7.38 


H0707 


0.091 


3.97 x 10 11 


7.44 


PG1211 


0.109 


2.21 x 10 11 


7.50 


PG1244 


0.191 


2.39 x 10 1() 


7.51 


PG1404 


0.081 


9.50 x 10 11 


7.62 
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4. Variability of NLSls 

The strength of the soft excess is correlated with 
the variability amplitudes; objects with strong soft ex- 
cesses show higher amplitude variability (Leighly 1999b). 
Hence, successful models for NLSls should account for 
the presence of soft X-ray variability. Leighly (1999a) 
summarize the ASCA observations of NLSls, finding sys- 
tematically large excess variance in NLSls. This fact 
can be, in principle, explained in terms of small black- 
hole masses, however, origin of giant-amplitude varia- 
tions observed in some NLSls, including IRAS 13224- 
3809 (Boiler et al. 1997, see also Otani et al. 1996) and 
PHL 1092 (Brandt et al. 1999), seems to require addi- 
tional explanation (see section 3.5). The lack of optical 
variability in IRAS 13224-3809 poses a problem (Young 
et at. 1999), although both of the soft and hard com- 
ponents (in the ASCA band) are known to vary nearly 
simultaneously in the same source (Otani et al. 1996). 
Relativistic boosting of otherwise modest variations may 
produce extreme variability of this source (Boiler et al. 
1997). 

The light curves of NLSls, especially those exhibiting 
extreme variations, have the following characteristics: 

(1) The light curves seem to be composed of numerous 
flares (or shots). 

(2) The flares do not have an identical time profile but 
show a variety of flare amplitudes and durations. 

(3) Temporal distribution of each flare is apparently ran- 
dom. 

Interestingly, these features are the same as those 
found in GBHC during the low (hard) and very high 
states (see, e.g., van der Klis 1995). Fluctuations in the 
soft X-ray component (of a disk origin) is observed only in 
the very high state, while the low-state disk has no soft 
component. In this sense, NLSls are more like GBHC 
during the very high state. 

To reproduce theses characteristics of the light curves, 
Kawaguchi et al. (2000) proposed that fractal magnetic- 
field structure spontaneously arises in an optically thin, 
advcction-dominatcd flow and produces sporadic mag- 
netic flare events, like solar flares, via magnetic recon- 
nection. In fact, dynamics of magnetic fields in ADAFs 
is known to be well represented by a cellular-automaton 
model based on the notion of the Self-Organized Critical- 
ity (Mineshige et al. 1994b, Kawaguchi et al. 1998). It is 
of great importance to note that both of optically thick 
and thin ADAFs may share common time-dependent 
properties, since they are both dynamical systems in the 
sense that radial accretion velocity, v r , is comparable to 
the free-fall velocity and sound speed, as long as a is not 
small. We, here, argue that the same MHD model can be 
relevant to the slim disk, thus explaining the variability 
of NLSls, as well. 

If fluctuations are of magnetic origin, large-amplitude 
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fluctuations indicate relatively large field energy com- 
pared to the radiation energy. Suppose that gas cloud 
of a unit mass falls onto a black hole. Since radiative 
cooling is efficient in the standard-type disk, energy that 
is emitted away is comparable to the release of the grav- 
itational energy, i? rac i <~ E glav (~ GM/r). As a result, 
the internal energy of gas should be much lower than the 
gravitational energy, -E gas <C -E grav . Since magnetic field 
energy will be, at most, the same order of the gas energy, 
^mag^^gas, we have £ mag < i?g rav ~ £ rad . If fluctu- 
ations occur due to sporadic release of magnetic energy 
via magnetic reconnection, this inequality explains small 
fluctuations in the standard disk. 

In the ADAF state, both optically thick and thin ones, 
radiative cooling is inefficient. Suppose again that gas 
cloud of a unit mass falls onto a black hole. A distinc- 
tion exists in what carries energy to a hole; it is internal 
energy of gas for the case of the optically thin ADAF, 
whereas it is trapped photons for the slim disks. In any 
cases, internal energy of gas or of trapped photons is 
of the order of gravitational energy which much exceeds 
radiation energy which is emitted away, (-Erad)trappcd > 
(-Erad)out- Since again £ mag < (£ rad ) tra p pod , moderately 
large magnetic energy compared with the energy radi- 
ated away is expected, -E m a g <; (-Erad)out- The larger m 
is, the more becomes the ratio of (E lad ) tlapped / (E lad ) out . 
Namely, large fluctuations in soft X-ray emission are in- 
evitable at high m, consistent with the observations of 
NLSls. 

In contrast, the standard-type disk (which may exist in 
BLSls) is stable and, hence, is not variable. Fluctuations 
from BLSls are thus due to variation of incident power- 
law radiation from disk coronae, which are presumably 
filled with magnetic fields and thus variable for the same 
reason mentioned above. In terms of variability, there- 
fore, the disk main body is passive in BLSls while it is 
active in NLSls. To elucidate the theory distinguishing 
these two, however, realistic modeling of the disk-corona 
structure is needed as future work. 

5. Outstanding issues 

5.1. Radiation transfer 

In a hot and dense layer with a very large Thomson 
depth (t), there are several important effects to be con- 
sidered in addition to Compton scattering effects. First, 
the photon diffusion time (~ Ht/c) may exceed the ac- 
cretion time. 

T dif ~ 4 x 10 4 S 4 (iJ/r)(r/10rs)M 5 s (10) 

and 

r acc ~ 5 x 10 2 (« r /0.1« ff )- 1 (r/10r s ) 3 / 2 M 5 s, (11) 

where £4 = E/10 4 g cm -2 and v r is the accretion (ra- 
dial) velocity. Then, the amount of photons which can 
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Table 3. Basic disk models and their properties. 
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disk models 



standard disk (r > 1) 



ADAF (r < 1) 

Qvis Qadv Qr&d 

H<r 
<10 10 K 

L oc (Mf 
power-law 

internal energy of gas 

^mag ~ -Egrav EvaA 

large 



ADAF (r > 1) 



energy balance 
disk height 

electron temperature 

ion temperature 

luminosity 

spectrum 

final energy form 

magnetic energy 

fluctuation 



Qvis Qrad ^ > Qadv 

H <r 

~ 10 6 M^(r/r s )-iK 

~ f0 6 M r 7'(r/r s )-3K 
L oc M 
blackbody 
emitted radiation 

^mag Lg ra v ~ L^rad 

small 



Qvis = Qadv > Q rad 

H<r 

~2x 10 6 M 5 _i (r/r s )"5K 

~2x f0 6 M 5 _ '(r/r s )"3K 
L oc In M 
blackbody 
trapped radiation 

Vag ^ Lg ra v 

large 



be emitted to space is limited so that the disk luminosity 
gets even lower, say below Le, even for large to. Such 
effects as those arising due to a finite diffusion time are 
omitted in the present study but this treatment is rea- 
sonable if efficient convective mixing of gas and radiation 
occurs so that photons generated on the equatorial plane 
can reach the disk surface on dynamical timescales. Such 
effects should be examined in future, however. 

Second, thermalization timescales may exceed the ac- 
cretion time for a <^ 0.03 (Beloborodov 1998). Then, disk 
temperature gets remarkably high, ~ 10 8 K even for the 
parameters of AGN, producing much higher-energy ra- 
diation. However, such high-energy emission spectra do 
not reconcile with the observed soft spectra of NLSls 
with soft excess of kT <~ lOOcV. Probably, again con- 
vective mixing of gas and radiation promotes rapid ther- 
malization between them. Or alternatively, the viscosity 
parameter could be small, a < 0.03, for which thermal- 
ization occurs within accretion time. 

Finally, radiation field is more likely to be anisotropic; 
i.e., radiation mostly goes out in the vertical direction to 
the disk plane, whereas matter approaches to the black 
hole along the disk plane. This much reduces flux to an 
observer viewed from the direction with large inclination 
angles. Sub-Eddington luminosity is expected even when 
the total L exceeds Le. 

5.2. Limit-cycle oscillations 

Since the standard and slim-disk branches form an S- 
shape in the thermal equilibrium curves (Abramowicz 
et al. 1988), there is a possibility of limit-cycle oscil- 
lations between the two branches (Honma et al. 1991; 
Szuszkiewicz, Miller 1998). Such oscillations, if occurs, 
will possibly explain quasi-periodic outburst behavior, 
which may explain transient AGN phenomena in oth- 
erwise normal (inactive) galactic nuclei (see, e.g., Ko- 
mossa, Bade 1999; Komossa, Greiner 1999), although re- 
peated outbursts will produce permanent AGN-typical 
NLR emission lines, which are not observed in some 



sources. Alternatively, the S-shaped curve may produce 
soft-hard transitions; i.e., transitions from the slim disk 
to the optically thin ADAF (Takeuchi, Mineshige 1998). 
These possibilities should be elucidated in future. 

5.3. Implications on the cosmological evolution 

It is often argued that for a BH to increase its mass 
by a factor of ~ e it takes is = M/Me ~ 10 8 yr, since 
the maximum accretion rate onto a BH is regulated by 
the Eddington rate, Me ~ 10Le/c 2 . Then, it takes quite 
a long time, ~ 2 x 10 9 yr, to make a 1O 9 M BH out 
of <~ IOMq BH. Much shorter time is enough, however, 
if supercritical accretion is allowed. Indeed, supercritical 
accretion is unavoidable to account for the present quasar 
population (e.g. Kauffmann, Haehnelt 1999). It might 
be that there used to be numerous NLSls in the past. 
The observational consequences need to be investigated. 

6. Conclusions 

To summarize, we discussed various aspects of the slim 
disk in the context of NLSls in comparison with the other 
disk models (see table 3). The followings are the main 
conclusions. 

1. The disk begins to exhibit the slim-disk properties, 
when the luminosity moderately exceeds Le- Then, 
substantial amount of radiation is expected to orig- 
inate from inside the marginally stable last circular 
orbit at 3rg. Therefore, a small rbb < 3r$ does not 
necessarily mean the presence of a Kerr hole. 

2. If we fit the spectrum of the slim disk with the disk 
blackbody model, the derived size of the X-ray emit- 
ting region will be r; n < 3r§ and the maximum tem- 
perature will be kT in ~ 0.2M 5 _1 ^ 4 keV. Further, p 
(= d In T e ff/d In r) decreases from 0.75 to <~ 0.5 as 
M increases. 
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3. All the ASCA data of NLSls fall on the region with 
supercritical accretion, m > 10, supporting that the 
disks in NLSls are likely to be slim disks. 

4. The slim disk can also produce large fluctuations 
because release of magnetic-field energy (which will 
produce fluctuations via magnetic flares) dominates 
over the energy release due to persistent emission. 
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